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Abstract 

grc4f is a Monte-Carlo package for generating e + e~ — * 4-fermion processes in the standard model. 
All of the 76 LEP-2 allowed fermionic final state processes evaluated at tree level are included in version 
1.1. grc4f addresses event simulation requirements at e + e~ colliders such as LEP and up-coming linear 
colliders. Most of the attractive aspects of grc4f come from its link to the GRACE system: a Feynman 
diagram automatic computation system. The GRACE system has been used to produce the computational 
code for all final states, giving a higher level of confidence in the calculation correctness. Based on the 
helicity amplitude calculation technique, all fermion masses can be kept finite and helicity information 
can be propagated down to the final state particles. The phase space integration of the matrix element 
gives the total and differential cross sections, then unweighted events are generated. Initial state radiation 
(ISR) corrections are implemented in two ways, one is based on the electron structure function formalism 
and the second uses the parton shower algorithm called QEDPS. The latter can also be applied for final 
state radiation (FSR) though the interference with the ISR is not yet taken into account. Parton shower 
and hadronization of the final quarks are performed through an interface to JETSET. Coulomb correc- 
tion between two intermediate W's, anomalous coupling as well as gluon contributions in the hadronic 
processes are also included. 



PROGRAM SUMMARY 

Title of program: grc4f vl.l 

Program obtainable from: CPC Program Library, 
Queen's University of Belfast, N. Ireland ( see ap- 
plication form in this issue ) and from ftp.kek.jp 
in directory kek/minami/grc4f. 
Computer for which the program is designed and 
others on which it is operable: HP9000 and most 
of the UNIX platforms with a FORTRAN77 com- 
piler 

Computer. HP9000; Installation: National Labo- 
ratory for High Energy Physics (KEK), Tsukuba, 
Ibaraki, Japan 

Operating system: UNIX Programming language 
used: FORTRAN77 

High speed storage required: 36 Mbyte Card image 
code: ASCII 

Key words: VF-boson, Z-boson, e + e~ colliders, 
four-fermion, ISR, FSR, QEDPS, Coulomb correc- 
tion, event generator, parton shower, hadroniza- 
tion, color base. 



Nature of physical problem 

Study of W boson physics and of background pro- 
cesses for new particle search at LEP-2 and be- 
yond. 

Method of solution 

The automatic amplitude generator GRACE is used 
to get the necessary helicity amplitudes for all 
the four-fermion processes. Specific corrections 
are then introduced to deal with the gauge boson 
width, the radiative corrections, the color assign- 
ment, the hadronization, the coulomb correction 
and the anomalous coupling. 
Typical running time 

The running time depends on the number of di- 
agrams of the selected process, on the required 
cross-section accuracy and on the applied cuts. 
For instance, the process e~i7 e ud takes 20 min- 
utes to reach a 0.5% accuracy on the total cross 
section integration with HP-735/99. 



1. Introduction 

The main purpose of the LEP-2 experiments is to measure the properties of the VF-boson 
with a high level of precision by its direct pair production. At the Z° peak (LEP-1 energy), 
two-fermion final states produced by the Z° decay are, by far, the most important processes. At 
higher energy, charged VU-boson being essensially created in pair at e + e~ colliders, the cross- 
section is dominated by four-fermion processes. However, other production mechanisms may 
also lead to four-fermion final states like Z° pair or two-photon like production. In addition, 
some four-fermion processes contribute heavily to the background for new particle searches due 
to their large missing energy. In this paper, a four-fermion event generator, grc4f , is presented. 
It provides a convenient way of computing cross-sections under the complex cuts, acceptance 
and resolution of the collider experiments. 

It is based on the GRACE [jl]] system which generates automatically the matrix element in 
terms of helicity amplitudes (supplied by the CHANEL library) for any process once the initial 
and the final states have been specified. At present, the system has been completed only at tree 
level in the framework of the standard model for electroweak and strong interactions. Fermion 
masses are non-zero and helicity information can be traced down to the final state particles. 
In addition, a kinematics library has been developed to cover the requirements of each process 
topology (weeding out multidimensional singularities) to get a faster convergence of the Monte- 
Carlo integration over the phase space. The grc4f package is actually a collection of 76 e + e~ — * 
4-fermion processes Q, consisting of 16 hadronic, 36 semi-hadronic and 24 leptonic processes, 
presented in a coherent and uniform environment. Flavor mixing is set to zero. The coupling 

^Jn version vl.l, final states with top quarks or those obtained by complex conjugation are not 
explicitly included. For generating i-quark final state processes, the user may try to change the c-quark 
mass to the top mass together with an interchange s «-» b, but the integration convergence is not warranted 
because t — > bW is a real process contrary to c — > sW. 
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between electron and Higgs boson is suppressed. Once a process has been selected, the total 
and differential cross sections are computed with the Monte-Carlo integration package BASES 
P]. Then SPRING [H], a general purpose event generator, provides unweighted events. 

The matrix element generated by GRACE corresponds to the genuine tree level process. How- 
ever, several corrections must be introduced to produce realistic cross sections. First of all the 
most basic correction, inherent to heavy boson production, is due to the finite width of the 
W or Z° boson to ensure a finite cross section. It should be noted that if one introduces the 
width in a naive way, it violates the gauge invariance and in some case may lead to a divergent 
cross-section. The e + e~ — > e~v e qq' cross section, for example, blows up at 9 e « 0. To cure this 
problem, the terms which satisfy the Ward identity are properly subtracted from the electron 
current which is connected to the t-channel photon Q. The computation of the constant and 
the running decay widths of and Z° is presented in the section 2.1. 

Two additional effects, although of different origins, can be, in first approximation, intro- 
duced into any selected process; radiative corrections and hadronization of the final quarks. 

For the radiative corrections two techniques are proposed in the program. The first one 
uses the well-known analytic form of the e ± structure function || and the second is based on 
QEDPS ||, a radiative correction generator producing an indefinite numbers of photons according 
to the parton shower algorithm in the leading-logarithmic (LL) approximation. Originally this 
algorithm has been developed to simulate QCD parton shower. One important point here is that 
QEDPS reproduces the radiative photon transverse momentum distributions. The details of these 
two methods shall be given in section 2.2. Final state radiation (FSR) for the electron and muon 
can also be generated with QEDPS. However the contribution of the interference between ISR and 
FSR is not yet included. Hence one should not use the FSR mode when the final state contains 
electron for kinematical regions where the t-channel photon exchange is the dominant diagram. 
The FSR for the quarks is not allowed in grc4f , because the QCD evolution scale is much shorter 
than the QED one, therefore the photonic correction on quarks might be meaningless. We 
recommend to use the FSR mode only for the case where the VF-pair or Z-pair (Zj) production 
diagrams are dominant. 

We assume that the hadronization of partons can be separated from the hard interaction 
studied here. Under this assumption, the calculation of cross sections is exact in grc4f . Final 
state parton hadronization is performed, in grc4f , through the LUND mechanism as implemented 
in JETSET|7j. A set of color bases is defined to select final state color flow. Possible ambiguities 
in this procedure as well as other QCD related issues are discussed in details in section 2.3. 

In order not to waste cpu time and disk space in the generation of events that will be 
eventually rejected by detector acceptance or threshold, a set of general cuts whose values can 
be set by the user have been implemented directly in the kinematics as presented in section 2.4. 

The Coulomb correction^] due to the exchange of a virtual photon between two intermediate 
W is important for non-relativistic W, namely close to the production threshold. It may reach 
4% of the total cross-section. The introduction of this effect is described in section 2.5. 

Anomalous coupling in the interaction among vector bosons is introduced in section 2.6. 

The basic input parameters used in the program are listed hereafter. Fermion masses are 
taken from the report of Particle Data Group ||, but can be set to zero if necessary. 

M w = 80.23 GeV, 
M z = 91.1888 GeV, 

T w = 3G f M^/(V8tt), G f = 1.16639 • 1(T 5 GeV~ 2 , 
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T z = 2.497 GeV, 
a(0) = 1/137.036, a(M w ) = 1/128.07, 
a s = 0.12. 

Here a(Q) is the QED fine structure constant at the energy scale Q and as is the QCD coupling. 
These values can be changed through control data. The weak mixing angle is calculated by 
sin 2 9 W = 1-M^/M§. 

The matrix elements are given in term of helicity amplitudes, it is therefore possible to select 
any helicity state configuration (initial and final). 

This paper is organized as follows. Theoretical details for the problems mentioned in the 
introduction will be given in section 2. The structure of the program is described in section 3. 
In section 4, all details about running the program are presented. Some improvement for future 
releases are summarized in the section 5. Three appendices present the parameters and options 
which can be changed by the user, the list of all processes and the program installation procedure 
followed by a run example. 

2. Theoretical aspects 

This section covers the treatment of the boson width, the introduction of the initial and 
final state radiation, the QCD related issues such as hadronization and the color problems, the 
description of kinematics and cuts, the Coulomb correction and the anomalous coupling of heavy 
vector bosons. 

2.1 Gauge boson width 

The 0{a) self-energy of the gauge boson, S^^(A: 2 ), generally satisfies 

Re£ (1) (M 2 ) = 0, 
ImS (1) (M 2 ) = -AiT(M 2 ), 

where M is the mass of the gauge boson and T(M 2 ) is the lowest order decay width. This 
induces the replacement of the gauge boson propagator by the following simple Breit-Wigner 
form: 

1 1 



k 2 — M 2 k 2 -M 2 + iMT(M 2 ) 

This form is used for the fixed width scheme in grc4f whenever the propagator appears with 
positive momentum squared. 

Another way is to take into account the energy dependence of the self-energy. Based on the 
observation that the contribution of fermion pairs to ImEW is k times a constant in any gauge 
as long as all the light fermion masses are neglected, the following energy dependent width has 
been proposed [1C]. 

k 2 

MT(M 2 ) -» -ImS«(fc 2 )| fermion pairs ^r(M 2 ). 

In grc4f this form is used as the default, but the constant width can be selected as well. 

As mentioned in the previous section, the introduction of the boson width may raise dif- 
ficulties, particularly when the electron is scattered in the very forward direction. Such final 
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states usually come from the generation of the so-called single W or Z process. To get rid of 
the divergent cross section, we apply the following method. 

With an arbitrary gauge parameter £ the electron current can be written as: 

l» = Hp')7 u u(p)lg^ + (Z-l)k^k u /k 2 }, (l) 
= u(p') llx u(p), (2) 

where p^{p'^) is the four- momentum of the initial (final) electron and kp, = p^ — p'^, the momen- 
tum of the virtual photon. After squaring the amplitude and averaging over spin states, the 
matrix element can be written as 



M 



J-'pu 1 > 



with 



L 



spin 



PtiP'u + PvPn + — 9p,v 



Let's assume that T^ u is gauge invariant. Then one can replace L^ u by 



(3) 
(4) 

(5) 



In this form the first term is responsible for the blow-up of the cross section. Thanks to the 
gauge invariance, one can further replace the vector p^ by: 



Pp. P/x = Pp. ~ (poAo) V 
By substituting P^ in Eq.(5) and dropping kn, one gets the final form 



L' 



T" 



AP^ + k'g^. 



(6) 



(7) 



It is known that a product of P with an arbitrary vector A, P • A, can be expressed by a sum of 
terms proportional to either m 2 , 1 — cos# e or sin# e . The part is expressed in terms of some 
momenta which are to be contracted with P^ or P u . Hence in the region 9 e ks 0, any product 
behaves effectively like k 2 , because 1 — cos(9 e vanishes almost like k 2 . If the current L"^ is used 



0. 



instead of the original one L^ v , the cross-section remains finite down to t 

2.2 Radiative corrections 

In the first approach, the simple electron structure function is used. The electron structure 
function at 0(a 2 ) || which is to be convoluted to the cross section for a primary process is 
given by 



D(x,s) 



1 + ^ + 

P 



/_9__C(2) 
V128 8 



2 V 







2 

4 V ' 32 
(2a/n)(ln(s/m 2 e ) 



4(1 + x) ln(l - x) + 
-1), 



x/3/2-1 



1 + 3x 2 
1 — x 



\nx + (5 + x) 



(8) 
(9) 



where s is the square of the total energy of the system, x is the momentum fraction of the 
electron, and C( n ) is Rieman's (^-function. It should be noted that, when compared with the 
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exact 0(a) calculation, the radiative corrected cross section obtained by this function does not 
contain the overall multiplying K- factor: 



q / 7T 2 1\ , . 

1 + - —-- = 1.006480 (10) 



7T 



In grc4f this factor is missing both for the structure function mode and QEDPS mode. If 
necessary, the final result can be multiplied by this factor for better accuracy. 

The QED Parton Shower approach, QEDPS, is primarily based on the fact that D(x,Q 2 ) 
obeys the Altarelli-Parisi equation: 

dD(x,Q 2 ) a f 1 dy 2 

-P + (x/y)D(y,Q ), (11) 



dlnQ 2 2tt J x y 



in the leading-log(LL) approximation [11 1. This is equivalent to the following integral equation: 

n rQ 2 dK 2 r l ~ e dti 

D(x,Q 2 ) = U(Q 2 ,Q 2 )D(x,Q 2 s ) + — / — Y U(Q 2 ,K 2 ) / ^-P(y)D(x/y, K 2 ), (12) 

where the small quantity e will be specified later. In these equations P{x) is the split function 
noted P+(x) when regulariz at x = 1. Q 2 is the initial value of Q 2 . For simplicity the fine 
structure constant a is assumed not running with Q 2 . The Sudakov factor II is given by: 

Q 2 dK 2 r 1 -* 



Q") = exp j Q2 J q dxP{x) j . (13) 

This is the probability that an electron evolves from Q' 2 to Q 2 without emitting hard photon. 
In other words, II already contains the soft photon component, which causes the change in the 
electron virtuality, and the loop correction contribution at all orders of perturbation. 

The integral form Eq.(|l2|) can be solved by iteration. It is clear that the emission of n 
photons corresponds to n iterations. Hence it is possible to regard the process as a stochastic 
mechanism suggesting the following shower algorithm E 



(a) Set Xf) = 1. The variable x^, is the fraction of the light-cone momentum of the virtual 

electron that annihilates. 

(b) Choose a random number r\. If it is smaller than U(Q 2 ,Q 2 ), then the evolution stops. If 

not, one can find the virtuality K 2 that satisfies rj = U(K 2 ,Q 2 ) at which a branching 
takes place. 

(c) Fix x according to the probability P(x) between and 1 — e. Then x\, is replaced by x^x. 

Go back to (b) by substituting K 2 into Q 2 and repeat until the evolution stops. 

Once a radiative process is fixed by this algorithm, each branching of a photon is a real process, 
that is, an electron with x, K 2 decays like: 

e~(x, —K 2 ) - e~(xy, -K' 2 ) + 7 (x(l - y), Q 2 ). (14) 

Here we have introduced a cutoff Qq to avoid the infrared divergence. The momentum con- 
servation at the branching gives: —K 2 = —K' 2 /y + Qq/(1 — y) + k|./(y(l — y)) which in turn 
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determines the photon transverse momentum relative to the parent, k 2 ^, from y,K 2 ,K' . This 
technique gives the distribution as well as the shape of x. 

The kinematical boundary y{K 2 + Qq/(1 — y)) < K' 2 , equivalent to k^ > 0, fixes e as 
e = Qq/K' 2 since K 2 <C K' 2 is expected. In ref.|jl3| the important role played by this e is 
discussed in more details. 

The above description of the algorithm concerns the case where either e~ or e + radiates 
photons when the axial gauge vector is chosen along the momentum of the other electron, 
namely e + or e~ . In the program, however, we use the double cascade scheme to ensure the 
symmetry of the radiation between e + and e~ jL4j]. These two are mathematically equivalent in 
the LL approximation. 

The two parameters Q 2 and Qq are given as follows; 

Q 2 s = m 2 e e = m 2 e x 2.71828 • • • , Q 2 = 1(T 12 GeV 2 . (15) 

Q 2 is defined so as to include the constant term — 1 of /3 in such a way that f3 = (2a/7r)(ln(s/m 2 ) — 
1) = (2a/n) ln(s / (m 2 e)) . Since Qq is unphysical, any observable should not depend on it. It 
has been checked that increasing Qq up to 0(m 2 /10) leaves the result unchanged within the 
statistical error of the event generation ||. 

This scheme can be applied for the radiation from the final state charged particles as well. 
The lower bound of the virtuality integration, Q 2 , is now the mass square of the final particle 
instead of the initial electron mass. The upper bound is the four- momentum squared of the 
lepton pair. Here we assume that the lepton-pairs are created from the gauge boson. The FSR 
should be used for the processes in which the VF-pair or Z-pair (Zj) production diagrams are 
dominating, but not for the processes such as multi-peripheral two-photon like diagrams. 

2.3 QCD related issues 
2.3.1 Color bases 

We calculate the matrix element using GRACE which may be written as 

M = \J2Tj\ 2 . (16) 

Each amplitude Tj includes a color factor. The color indices of the external particle must be 
summed in the final state and averaged in the initial state. 

If a diagram has four gluon vertices, the diagram is separated into three pieces, the so-called 
s-, t- and u-channels, respectively, for each vertex. Let's now consider the Tj's not as amplitude 
but as a component of an amplitude after this decomposition. 

For each Tj, the color factor can be factorized out and expanded on a set of color bases, 




where Wj s are numbers. Then the matrix element is given by 

k,k' j 

Therefore, for the four-gluon vertices, the colored part of an amplitude consists of three- 
gluon vertices (—ifabc), quark-gluon vertices (t a ) and non-colored vertices. The ghost-gluon 
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vertex bears the same color structure as the three-gluon vertex and do not require a specific 
treatment. 

The choice of the color bases and the technique used to introduce color flow at the event 
generation level are not all equivalent. Let's describe one of the simplest approach. The following 
algorithm has been introduced in GRACE. 

1. We consider a process in which n q quark pairs and n g gluons exist among external particles. 
For any pair of amplitudes, Tj and t1, the color index of an external gluon is contracted 
between them. Making use of this, each gluon is converted into a pair of quark and 
antiquark creating a quark-gluon vertex. An overall factor 2 n Q is assigned to M from: 

5 ab = Tr t a t b /T R , TR = \- (19) 

2. External particles are n pairs of quarks and antiquarks where n = n q + n g . Quarks and 
antiquarks are denoted as 1, 2, 3, . . . , n and 1, 2, 3, . . . , n, respectively. The color base here 
has the form: 

Ck = Sij!^2j 2 ' ' ' $nj n (20) 

where ji,j2, ■ ■ ■ ,jn is a permutation of I, 2, 3, ... , n. The number of color bases is hence 
n\. 

3. Each three-gluon vertex is converted into a pair of quark loops by 

- ifabc = (-Tr t a t b t c + Tr t b t a t c )/T R . (21) 

4. Each gluon propagator is replaced by quark lines by use of the Fiertz transformation: 

(t a )ij5 ab (t b )ki = --^-Sij5ki + -SiiSjk, (Nc = 3) (22) 

5. After the above procedures, only quark lines remain. For each closed quark loop, a factor 
Tr 1 = 3 is assigned. 

When the calculation is done in the covariant gauge, the diagrams including external ghost 
particles, if any, are to be considered separately^. 

A different approach can be followed without gluon-quark conversion. For example, let 
us denote quarks, antiquarks, and gluons as 1, 2, 3, . . . , n g , 1, 2, 3, . . . , n q , and 1, 2, 3, . . . , n g , 
respectively. Instead of Eq.(20), the color bases can be represented by a product of the following 
objects, 

Sin, (^)yi, (^)iii, ••• 
where ji,j2, ■ ■ ■ and ki, k?,,... are permutation of 1, 2, 3, ... , n q and 1, 2, 3, ... , n g , respectively. 

2.3.2 Interface to hadronization 

As mentioned in the previous section, there are many possible choice for selecting a color 
base. As far as the total cross section is concerned they give the same result. However, the 
choice of the color base and the technique used to introduce color flow at the generation level 
are not equivalent and may induce noticeable discrepancies after hadronization. 



2 This does not happen in grc4f . 
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We take the string picture used in JETSET as a concrete example of hadronization. The color 
flow pattern can be assigned to each color base intuitively. However, an event is generated by 
matrix element which is the square of the amplitude, i.e., a linear combination of color bases. 
For interference terms, assignment of a color flow pattern is ambiguous. One possible solution 
is to introduce a set of orthogonal bases. 

To illustrate this point, we consider the four-quark final state, qi&q.stlA which is produced by 
grc4f . We label them 1, 2, 3, and 4. Among possible candidates, let us consider the following 
three types of color bases. 



1. Primitive base 



2. Orthogonal base 



3. Extended orthogonal base 



<U2#34 
<5l4<>32 



c{ o) =d + C 2 
C ( 2 o) =Ci-C 2 



— Cl — 5i234 



(23) 



(24) 



=C 2 -6 12U (25) 

= <5i234 

Here £1234 = 1 only when all four indices are equal and = otherwise. 

For each primitive base, one can assign a configuration of two strings naturally. They are, 
however, not diagonal; the interpretation of interference term is obscure. In the second case, the 
orthogonal bases, there is no interference term. However, each base do not relate directly 

to the string picture and each cjf* contains |Ci| 2 and | C 2 1 2 with equal weight. The same thing 
happens in elementary quantum mechanics: The information on left/right circular polarization 
is lost after the light goes through a linear polarization filter. In the third example, the first two 
can be assigned to two strings of different colors. The third one is a kind of 4-quark loop 
worm. They are orthogonal to each other but they do not fit well with the JETSET approach. 

The situation does not change for the case with more external particles, and for both cases 
with or without the gluon conversion. We have to say that there is some distance between the 
color bases and the definition of color singlet sub-system (e.g., strings) in the hadronization 
model. 

In the present scheme of grc4f , we use the primitive base case and the color flow is chosen 
for each event stochastically with weight CkC\\Fk\ 2 '. 

2.3.3 Options 

There is an option to include the overall QCD correction factor. Here the cross section is 
multiplied by a simple factor (1 + as/it) for each quark vertex. 

When the final state is a four-quark state, gluon exchange diagrams are included. It is 
possible to remove these by setting control data. 

2.4 Kinematics and cuts 

For the phase-space integration by the adaptive Monte Carlo method, a proper treatment 
of the singular behavior of the amplitude is necessary. It may come from: 
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1. s-channel W, Z resonances and 7 propagators. 



2. i-channel electron and electron neutrino propagators in Z-pair (Zy, 77) production and 
in TV-pair production. 

3. The electron (positron) forward scattering in two-photon processes. 

4. Mixed resonance of TV-pair and Z-pair (Z-y, 77) processes, such as udud process. 

5. Identical particles in the final state, such as /i + n~ ji + fi~ . 

grc4f includes the mapping routine from the eight integration variables to the four-momenta 
of final particles with proper treatment of these singular behavior. (A convergence of the nu- 
merical integration for the forward scattering of the two-photon processes is not so good.) Higgs 
boson diagrams are included in the amplitude, however the kinematics does not treat the Higgs 
boson resonance yet. (For cross section with very forward electrons or with Higgs boson res- 
onance, contact the authors 0.) For the identical particles in the final state, some momentum 
ordering is assumed in the kinematics routine, then, for instance, the first particle has always 
larger momentum than the second particle. 

One can apply some experimental cuts for the phase space integration; energy and angle cuts 
on all final particles and a invariant-mass cut on any pair of final particles. These parameters 
can be specified by the user. For comparison with other programs, grc4f has the canonical 



cut option used in ref.[15], which is; 

1. the energy of charged leptons must be greater than 1 GeV. 

2. the polar angle of charged leptons must be between 10 and 170 degree. 

3. the energy of quarks must be greater than 3 GeV. 

4. the opening angle between charged leptons and quarks must be greater than 5 degree. 

5. the invariant mass of quark pairs must be greater than 5 GeV. 

Here, the charged leptons include r's too. 

2.5 Coulomb correction 

This effect has been originally discussed for the off-shell TV-pair production || corresponding 
to the three resonant diagrams. The set of these diagrams cannot satisfy the gauge invariance 
and has no physical meaning in principle. In grc4f , however, all relevant diagrams are taken 
into account and the gauge invariance is restored from the contribution of the non-resonant 
diagrams (besides an effect from the finite width of IV boson). To maintain the invariance we 
adopt the following prescription: the Coulomb factor is applied to the minimum set of gauge 
invariant diagrams containing the W-pair production^L^. Hence the Coulomb multiplicative 
factor appears in some diagrams even without TV-pair. The following formula is used: 



air 

&Coul = O gauge ^'"■'i^ 



1 2 (\p M + A \ 2 -P 

1 arctan' 



7T V 2/?Im/3 



M 



3 E-mail address:grc4f Ominami . kek . jp 
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where 

P = - g ^s 2 - 2s(k% + kl) + (4 - kl ) 2 , 
Pm = y/l ~ 4M 2 /s, M 2 = Myy — iMwT\y, A = \ k +~ k -\ 

and — 7r/2 < arctany < tt/2. Here /3 is the average velocity of the W boson in its center-of-mass 
system. Two squared momentum k\ and k 2 _ are the virtualities of the intermediate W bosons. 

2.6 Anomalous coupling 

In the program, the anomalous coupling of heavy boson is available for the convenience of 
the user who may be interested to study such an effect, although there is no definite and reliable 
model beyond the standard model at present. The program includes only those terms which 
conserve C and P invariance which correspond to the following effective Lagrangian[|l6|] ; 

L eff = -igv{Wl v W»V v -W^W^-igvnvWlW^ - i gv ^_wl w»V Xv , 

m w 

= dpW v -d v Wn, v^ = d^v v -d u v^ v = z Q oi 1 , gv = 

Here Ky and Ay stand for the anomalous couplings parameters being 1 and 0, respectively, in 
the standard model. Deviation from these values corresponds to the introduction of anomalous 
coupling. 



-e V = 7 

-ecotflvK V = Z° 



3. Structure of the program 

The generator grc4f enables us to generate any e + e~ — > 4-fermions events with or without 
radiative corrections. This requires different calculation in some part of the program. In addition 
there are many options covering theoretical and experimental requirements. Since all program 
components are distributed as source code, users can set all options by editing the relevant 
subprograms directly. However, an interface program grc4f has been created to lighten the 
user's burden. It selects and/or corrects the components which are affected by the various 
options and create a "Makefile" according to the user requirements. This procedure is called 
the source generation phase. 

In the integration step, the matrix element of a selected process is integrated over the phase 
space by the subprogram BASES, which gives the effective total and differential cross sections and 
the probability distribution used in the event generation phase. There, the subprogram SPRING 
samples a point in the phase space and test if it can be accepted as a new event according to its 
probability. When an event is accepted the program control returns from SPRING to the main 
program, where further analysis is performed by using the resultant four-momenta of the final 
state particles. 

There are therefore, three steps in the generator grc4f , the first is the source generation, 
the second is the integration, and the third is the event generation. 

In addition to the user interface program grc4f , the following programs are available in the 
grc4f package: 

i) The main programs MAINBS and MAINSP, and all program components for the integration 
and event generation steps. 
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ii 



The interface programs to CERNLIB (GRC2CL) and to JETSET (GRC2SH). 



iii) 76 function programs FUMCs for processes e + e — > 4-fermions, each of which calculates the 
numerical value of the differential cross section for each process. 

iv) The kinematics subprograms KINMOQ and KINEMO for the processes with radiative correc- 
tion using the QED parton shower model or the electron structure function, respectively. 

v) The library CHANEL used to calculate the numerical values of basic components of Feyn- 
man diagrams in terms of helicity amplitudes. 

vi) The numerical integration and event generation program package BASES/SPRING v5 . 

The relationship among these program components and their function are presented in the next 
two sections. 

3.1 Source generation step 

The user interface program grc4f reads the parameters from control data prepared by users, 
which contains process selection, type of radiative corrections, physical options, experimental 
cuts, etc., and generates all necessary components: 

i) A main program MAINBS for the integration. 

ii) Interface programs GRC2CL and GRC2SH. 

iii) Four initialization subprograms USRPRM, MODMAS, KINIT, USERSP, and 

iv) a "Makefile". 



3.2 Integration step 

Before starting the numerical integration the main program MAINBS invokes an initialization 
subprogram USERIN, where the following subprograms are called in this order: 

USRPRM : To define the set of so-called "canonical cut" |l5| authorized by 

LEP200 working group and some additional optional parameters. 
SETMAS : To set masses and decay widths of particles. 
MODMAS : To alter the default values of all parameters defined in SETMAS. 
AMP ARM : To set the coupling constants and others parameters. 
KINIT : To set the parameters for the integration, kinematics, cuts, his- 
tograms etc. 

The subprograms SETMAS and AMPARM are generated by the GRACE system. There are no 
consistency checks among the constants, e.g. Mz, M\y and sin 9\y, so all modifications on these 
parameters in the subprograms USRPRM, MODMAS and KINIT are under the user responsibility 
alone. 

The integration program BASES calculates the scattering cross section by sampling the func- 
tion FUNC on the allowed phase space segmented by an self adapted grid where finer cells are 
clustered on the high gradient zones. This is an iterative process running until either the maxi- 
mum number of allowed iteration is reached or the required accuracy is obtained. In the function 
FUNC, the kinematics subroutine KINEMO or KINMOQ is used to map the integral variables with 
the four-momenta of the final state particles. KINEMO is used for reaction with no radiative 
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corrections or those involving initial radiation treated with the structure function techniques. 
KINMDQ is called for processes where radiative corrections are computed with the QED parton 
shower method. 

It is recommended to check the integration result carefully, especially the con- 
vergency behaviors both for the grid optimization and integration steps. When the 
accuracy of each iteration fluctuates, iteration by iteration, and, in some case, it may jump up 
suddenly to a large value compared to the other iterations, the resultant estimate of integral 
may not be reliable. There are two possible origins of this behavior; too few sampling points 
or/and an unsuitable choice of the kinematics. 

After the numerical integration by BASES, the subprograms BSINFO and BHPLOT are called 
to print the result of integration and the histograms, respectively. Before terminating the inte- 
gration procedure the probability distribution obtained by the integration can be saved in a file 
by invoking BSWRIT, which is then used for the event generation by SPRING. 

3.3 Event generation step 

After integrating the differential cross section and saving the probability distribution, the 
main program, MAINSP, handles the event generation program. The subprogram BSREAD is 
invoked to restore the probability distribution and then the subprogram USERIN is called. Each 
call to SPRING generates one event by sampling a point in the phase volume. It calculates first 
the differential cross section at that point using the same function FUNC seen in the integration 
phase and returns the weight of this sampling point. A weight-one event is finally generated 
using the usual unweighting technique. When an event is generated, SPRING returns the particle 
types and four momenta. The event information is stored in the labelled common LUJETS by 
calling the subprogram SP2LND. Here, the information for the color connection to be referred by 
JETSET is also supplied. Among the final states in Appendix B, those in Table 3 have non-trivial 
color flow which is determined as is described in the section 2.3. 

Then parton shower and hadronization of quarks and gluons can be performed by calling 
LUSHOW and LUEXEC. At the end of the event generation, the routine SPINFO and SHPLOT are 
invoked successively for printing event generation information and histograms. 

4. How to run the program 

The user should first prepare the control data to define the process, the options, flags and 
the experimental cuts. The user interface program grc4f takes this control data (let's call it 
control .data) as an input. 



process 


= eNEuD 


energy 


= 190. OdO 


canon = 


= yes 


type 


= tree 


massive 


= yes 


coulomb 


= no 


anomal 


= no 


qcdcr = 


= no 


end 





The first line specifies the process to be calculated and the second is the center of mass energy 
in GeV unit. The others are options, whose meanings are given in Appendix A. Then the user 
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may type: 



% grc4f < control . data 

If the message "syntax error" is returned, no files will be generated and the contents of the control 
data file must carefully be checked. After a successful completion, the following messages should 
be returned: 

Process is "eNEuD" 
Energy is "190.0d0" 
CANON <yes> 
MASSIV <yes> 
CDULMB <no> 
ANDMAL <no> 
QCDCR <no> 
bye-bye 

absolute directory name is /home/grc4f /prc/elNEuqDQ 

cd /home/grc4f /prc/elNEuqDQ 

make integ 

integ 

make spring 
spring 



According to the parameters given in control data, the files, i.e. usrprm.f, modmas.f, kinit.f, 
mainbs.f, grc2cl.f 1 grc2sh.f, usersp.f and Makefile, are generated in a specified subdirectory 
( elNEuqDQ in this case). 

According to the last four lines in the message, users can proceed the calculations as follows: 

i) Change directory by 

7, cd /home/grc4f/prc/elNEuqUD 

ii) Create an executable integ for the integration by typing: 
7o make integ 

iii) Numerical integration is actually performed by typing: 
To integ 

The results of integration step are displayed as well as written in an output file 
bases. result. The total cross section in pb with the error are displayed at the last 
row, under Cumulative Result, in the table of the Convergence Behavior for the 
Integration step. The differential cross sections as a function of the energy and angle 
of each particle and invariant masses of any two final particles will also be printed. The 
probability distribution is written in a file bases. data which will be used in the event 
generation step by spring. 
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) Before running the event generation, users should edit usersp . f to set additional param- 
eters like the number of events and call user's own analysis routines. 

The following is the structure of mainsp . f , where four- momenta of all particles are stored 
in the common/lujets/ in the JETSET format by calling subprogram sp21nd in the event- 
loop: 

implicit real*8(a-h,o-z) 



real*4 p,v 

common/lujets/n,k(4000,5) ,p(4000,5) ,v(4000,5) 



call usersp (mxtry ,mxevnt) 



do 100 nevnt = 1, mxevnt 

call spring ( func, mxtry ) 



* 

call sp21nd 

* 

* 

* ===================== 

* ( user_analysis ) 

* ===================== 

* 

100 continue 



stop 
end 

) Create an executable spring for event generation by typing: 

°/ make spring 
) Start the event generation by typing: 

°/ spring 

Information on the event generation will be written in the spring. result file. Users 
should pay special attention to the histograms. The distributions of generated events are 
superimposed by "0" on the original histograms by BASES. These two distributions should 
be consistent within the statistical error of the generation. For the detail of the output 
files of BASES and SPRING, user can consult the Ref.||. 
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5. Summary 

The generator grc4f enables us to calculate the effective cross section and to generate events 
for one of 76 e + e~ — * 4-fermions processes listed in appendix B. It is dedicated to the LEP and 
future linear collider physics studies. The produced quarks can be hadronized according to 
JETSET. Also processes with initial and final radiations can be generated in terms of the electron 
structure function or the QED parton shower methods, though the inclusion of the interference 
between the initial and final radiations requires further study. 

There still remain important problems to be solved and further necessary improvements. 
Among them it will be desirable to extend the program so as to produce several final state 
processes in a single run like e + e~ — > 4-quarks. In order to get more precise QED corrections 
an implementation of a complete next-to-leading logarithmic calculation is unavoidable. 

In this version, we have assumed that the hadronization of partons takes place independent of 
the hard interaction which produces partons. However a very important contribution arises from 
diagram where a gluon is exchanged between quarks produced in the two W decays. Taking 
into account this effect implies higher order calculations or the implementation of a specific 
phenomenologic correction [[l5|]. 

Finally the present version provides only 76 processes, the missing processes will be prepared 
soon by using GRACE. 
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Appendix A. Option table in Control data 

In the table below, the default values are underlined and the relation between commands and 
variable/array in Fortran sources is also described. Variable names are written in bold letters 
and filenames in italic. 

i) Process selection. 

Process = eNEuD 

abbreviation of process name 

This specifies the subdirectory name, where process related subroutines are stored. Tables 
1, 2, and 3 in Appendix B shows the abbreviation of process names and the subdirectory 
names where they are stored are listed. 



ii) 



Center of mass 


energy. 


energy 


190. dO 




CMS energy in GeV 




w in kinit.f 


Global options 


helicityl = 


average, left, right 




Helicity state for the initial electron. 


helicity2 = 


average, left, right 




Helicity state for the initial positron. 


helicity3 = 


sum, left, right 




Helicity state for 3rd particle. 


helicity4 = 


sum, left, right 




Helicity state for 4th particle. 


helicity5 = 


sum, left, right 




Helicity state for 5th particle. 


helicity6 = 


sum, left, right 




Helicity state for 6th particle. 



type 


tree, sf, qedpsi, qedpsif 




Type of calculation: 




Without radiation(tree), ISR with structure func- 




tion^), ISR with QEDPS(qedpsi) and ISR and FSR 




with QEDPS (qedpsif). 




jqedps = 0, isr = without radiation(tree). 




jqedps = 0, isr = 1 for sf. 




jqedps = 1, ips = 1 for qedpsi. 




jqedps = 1, ips = 2 for qedpsif. 




jqedps in usrprm.f, isr and ips in kinit.f. 
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canon 


yes, no 




Apply canonical cuts or not: 




j canon = 1 or in usrprm.f. 


massiv 


yes, no 




Quarks are massive or massless: 




jqmass = 1 or in usrprm.f. 


width 


run, fixed 




Running width or fixed width: 




j width = or 1 in usrprm.f. 



) Physical options 



coulomb 


yes, no 




Coulomb correction: 




jcolmb = 1 or in modmas.f. 


qcdcr 


yes, no 




Include overall QCD correction factor: 




j qcdcr = 1 or in modmas.f. 


gluon 


yes, no 




Include diagrams with gluon exchange: 




j gluon = 1 or in modmas.f. 



anomal = yes, no 

jano3v = 1 or in modmas.f. 
Anomalous coupling: 
l.OdO 

O.OdO 

l.OdO 

K Z 0. 

O.OdO 

X z o. 

See section 2.6 for the definition of above variables, 
if "anomalous = no" is specified, these options affect 
nothing. 

) Masses, widths and a (only if "canon = no"). 



ankaa 
anlma 
ankaz 
anlmz 
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amw 


= 80.23D0 

W mass in GeV: amw in modmas.f. 


agw 


= 2.033D0 

W width in GeV: agw in modmas.f. 


amz 


= 91.1888D0 

Z mass in GeV: amz in modmas.f. 


agz 


= 2.4974D0 

Z width in GeV: agz in modmas.f. 


alphai 


= 128.07D0 

a~ 1 :alphai in modmas.f. 


alpha_s 


= 0.12D0 

a s :alpha_s in modmas.f. 



vi) Experimental cuts (only if "canon = no"). 

The numbering convention of particles follows the GRACE scheme, where the initial electron 

and positron are 1st and 2nd, respectively, and the four final particles are numbered 3, 4, 

5, 6. In the process name of Table 1, 2, and 3, the order of particles corresponds to this 

numbering convention. For instance, in the process, e~u e ud, the 3rd is e~, the 4th is i7 e , 

the 5th is u and the 6th is d. 
_____ _ __________ 

Angle cut for 3rd particle in degree (backward- 
angle, forward-angle) . 

coscut (1 : 2, l)=cos (thecut3) in kinit.f. 
thecut4 = 180.d0 , 0.0d0 

Angle cut for 4th particle in degree (backward- 
angle, forward-angle) . 

coscut (1 : 2,2) =cos (thecut4) in kinit.f 
thecut5 = 180.d0 , 0.0d0 

Angle cut for 5th particle in degree (backward- 
angle, forward-angle) . 

coscut (1 : 2,3) =cos (thecut5) in kinit.f. 
thecut6 = 180.d0 , 0.0d0 

Angle cut for 6th particle in degree (backward- 
angle, forward-angle) . 

coscut (1 : 2,4) =cos (thecut6) in kinit.f. 

Instead of giving a numerical value, the user can use the keywords as below: amass 1(0 
has the mass for i-th particle and w is the CM energy. 
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vii) 



engcut3 


amassl (3) ,w 

Energy cut for 3rd particle (min.,max.) 
engyct (1 : 2 , 1) in kinit.f 


engcut4 


amassl (4) ,w 

Energy cut for 4th particle (min.,max.) 
engyct (1 : 2,2) in kinit.f 


engcut5 


amassl (5) ,w 

Energy cut for 5th particle (min.,max.) 
engyct (1 : 2 ,3) in kinit.f 


engcut6 


amassl (6) ,w 

Energy cut for 6th particle (min.,max.) 
engyct (1 : 2,4) in kinit.f 




ivmcut34 


amassl (3)+amassl (4), w-(amassl(5)+amassl(6)) 

Invariant mass cut for 3-4 particles(min.,max.) 
amasct (1 : 2 , 1) in kinit.f 


ivmcut56 


amassl (5)+amassl (6), w-(amassl(3)+amassl(4)) 

Invariant mass cut for 5-6 particles(min.,max.) 
amasct (1 : 2 ,2) in kinit.f 


ivmcut35 


amassl (3)+amassl (5), w-(amassl(4)+amassl(6)) 

Invariant mass cut for 3-5 particles(min.,max.) 
amasct (1 : 2 ,3) in kinit.f 




ivmcut46 


amassl (4)+amassl (6), w-(amassl(3)+amassl(5)) 

Invariant mass cut for 4-6 particles(min.,max.) 
amasct (1 : 2 ,4) in kinit.f 


ivmcut36 


amassl (3)+amassl (6), w-(amassl(4)+amassl(5)) 

Invariant mass cut for 3-6 particles(min.,max.) 
amasct (1 : 2 ,5) in kinit.f 


ivmcut45 


amassl (4)+amassl (5), w-(amassl(3)+amassl(6)) 

Invariant mass cut for 4-5 particles(min.,max.) 
amasct (1 : 2 ,6) in kinit.f 


Parameters for integration. 


itmx 


= L 15 

Iteration numbers: itmxl, itmx2 in kinit.f 


acc 


= 0.1, 0.05 

Accuracies in %:accl, acc2 in kinit.f 


ncall 


= 40000 

Sampling points: ncall in kinit.f 



viii) Parameters for event generation. 
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50 

Maximum trial numbers: mxtry in usersp.f 
10000 

Maximum event numbers: mxevnt in usersp.f 
yes, no 

Hadronization with JETSET0 for event generation: 
If yes, then a statement, call luexec is added in 
grc2sh.f. In this case, program JETSET is necessary, and 
one should specify JETSET file location in this control data, 
jetset = /home/jetset/jetset74.o 

Object file name for JETSET, if hadron = yes 
qccLshower = yes, no 

QCD parton shower with JETSET|?J for event generation: 
If yes, then a statement, call lushow is added in 
grc2sh.f. In this case, program JETSET is necessary, and 
one should specify JETSET file location in this control data. 





Those parameters affect Makefile. 




ix) HBOOK interface. 






cernlib = yes, no 

If yes, then the bases. hbook in integration and 
spring. hbook in event generation will be generated. 




Those parameters affect mainbs.f, grcBcl.f and Makefile. 




x) End of description. 






end 






After the command end any command is neglected. 





Appendix B. Process table 

Here, the 76 processes included in grc4f are listed. In the heading, 'abbrev.' and 'dir.' stand 
for the abbreviated name used in the control card and directory name where the generated code 
is stored, respectively. 



mxtry 

mxevnt 

hadron 
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process 


abbrev. 


dir. 


process 


abbrev. 


dir. 


+ - - 

6 J ^ 5 ^ ? 6 


eNEEne 


elELneNE 


e~ , Vv. u~*~ . v., 


eNEnmMU 


elNEMUnm 


e — . ZA,. T^. ZA- 


eNEntTAU 


elNETAnt 


U,, . . U~ . 9,, 

//3 3 3 /X 


muMUnmNM 


muMUnmNM 


Z/r- . . T. ZA- 


tauTAUntNT 


taTAntNT 


3 it 3 ' 3 T 


muNMTAUnt 


muNMTAnt 




eEeE 


elELelEL 


^ 3 ^ 3 3 A* 


eEmuMU 


elELmuMU 




eEtauTAU 


elELtaTA 


) r ! r )r 


muMUmuMU 


muMUmuMU 


r, r+,r, r+ 


tauTAUtauTAU 


taTAtaTA 


Zi _ . U . T. T 

p. 3 ^ 3 ' 3 ' 


muMUtauTAU 


muMUtaTA 




eEnmNM 


elELnmNM 




eEntNT 


elELntNT 




neNEmuMU 


neNEmuMU 


Vp.Vp, T, T + 


neNEtauTAU 


neNEtaTA 


V T , 9 T ,{1~,{1+ 


ntNTmuMU 


ntNTmuMU 




nmNMtauTAU 


nmNMtaTA 




neNEneNE 


neNEneNE 


v e ,9 e ,v^,9^ 


neNEnmNM 


neNEnmNM 


V e ,9 e ,V T ,9 T 


neNEntNT 


neNEntNT 




nmNMnmNM 


nmNMnmNM 




ntNTntNT 


ntNTntNT 


V T ,9 T ,V^,9^ 


nmNMntNT 


ntNTnmNM 



Table 1 Leptonic processes in grc4f . 



process 


abbrev. 


dir. 


process 


abbrev. 


dir. 


e~, 9 e , u, d 


eNEuD 


elNEuqDQ 


^ 3 3 5 


eNEcS 


elNEcqSQ 


H~ ,9^,u,d 


muNMuD 


muNMuqDQ 


ji~,9^,c,s 


muNMcS 


muNMcqSQ 


r, 9 T , u, d 


tauNTuD 


taNTuqDQ 


t,9 t ,c,s — bar 


tauNTcS 


taNTcqSQ 


e~ ,e + ,u, u 


eEuU 


elELuqUQ 




eEcC 


elELcqCQ 


e~ , e + , d, d 


eEdD 


elELdqDQ 


C ^ 6 S j S 


eEsS 


elELsqSQ 


e~ , e + , b, b 


eEbB 


elELbqBQ 


H~ , [i + ,u, u 


muMUuU 


muMUuqUQ 




muMUcC 


muMUcqCQ 


T, T + ,U,U 


tauTAUuU 


taTAuqUQ 


r, t + ,c,c 


tauTAUcC 


taTAcqCQ 


fx" , fi + , d, d 


muMUdD 


muMUdqDQ 




muMUsS 


muMUsqSQ 


fi~,fi + ,b,b 


muMUbB 


muMUbqBQ 


r, r + ,d,d 


tauTAUdD 


taTAdqDQ 


T, T + ,S,S 


tauTAUsS 


taTAsqSQ 


t, r + ,b,b 


tauTAUbB 


taTAbqBQ 


v e ,9 e ,u,u 


neNEuU 


neNEuqUQ 


v e , 9 e , c, c 


neNEcC 


neNEcqCQ 


v e , 9 e , d, d 


neNEdD 


neNEdqDQ 


3 3 ^ 3 ^ 


neNEsS 


neNEsqSQ 


u e ,9 e ,b,b 


neNEbB 


neNEbqBQ 


V^,9^,U,U 


nmNMuU 


nmNMuqUQ 


V/j,, IS/j,, c, c 


nmNMcC 


nmNMcqCQ 


v T , 9 T ,u, u 


ntNTuU 


ntNTuqUQ 


5 Z-^ -j- . C 


ntNTcC 


ntNTcqCQ 


v^,9^,d,d 


nmNMdD 


nmNMdqDQ 


^/i3 ^/i3 ^3 ^ 


nmNMsS 


nmNMsqSQ 


v^,9^,b,b 


nmNMbB 


nmNMbqBQ 




ntNTdD 


ntNTdqDQ 




ntNTsS 


ntNTsqSQ 


u T , 9 T , b, b 


ntNTbB 


ntNTbqBQ 



Table 2 Semi-hadronic processes in grc4f . 
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process 


abbrev. 


dir. 


process 


abbrev. 


dir. 


u, d, d, u 


uDUd 


uqDQdqUQ 


esse 


cCSs 


cqSQsqCQ 


u, d, s,c 


uDsC 


uqDQsqCQ 


U, u, u 


uUuU 


uqUQuqUQ 


c, c, c, c 


cCcC 


cqCQcqCQ 


d, d, d, d 


dDdD 


dqDQdqDQ 




sSsS 


sqSQsqSQ 


b, b, b, b 


bBbB 


bqBQbqBQ 


u, u, c, c 


uUcC 


uqUQcqCQ 


u, u, s, s 


uUsS 


uqUQsqSQ 


u, u, b, b 


uUbB 


uqUQbqBQ 


c, c, d, d 


cCdD 


cqCQdqDQ 


c, c, b, b 


cCbB 


cqCQbqBQ 


d, d, s, s 


dDsS 


dqDQsqSQ 


d, d, b, b 


dDbB 


dqDQbqBQ 


s, s, b, b 


sSbB 


sqSQbqBQ 



Table 3 Hadronic processes in grc4f . 



Appendix C. Installation 

The source code is available by anonymous ftp from ftp.kek.jp in the directory 
kek/minami/grc4f . The grc4f package contains the complete set of Fortran sources for 76 
processes, the three libraries, i.e., BASES/SPRING, CHANEL and utilities for kinematics. Those 
source codes are written in FORTRAN77. In addition, grc4f provides the interface program 
to generate a few Fortran source files according to the control data specified by the user. This 
program is written in C, YACC and LEX. grc4f has been developed on HP-UX, but should run 
on any UNIX platform with a fortran complier. 

The procedure of installation is as follows: 

1. Editing Makefile. 

The following macros in Makefile should be taken care of by users themselves. For example, 
in the right hand side of GRC4FDIR the directory name where grc4f is installed should be 
given, and for FC and FOPT the relevant compiler name and option for your system should 
be given. The other macros can be left as they are. 



GRC4FDIR = 


directory where grc4f are installed. 


PRCDIR 


directory where process files are installed. 




(default is $(GRC4FDIR)/prc.) 


LIBDIR 


directory where libraries are installed. 




(default is $(GRC4FDIR)/lib.) 


BINDIR 


directory where an executable is installed. 




(default is $ (GRC4FDIR) /bin.) 


MACHINE 


[hpux | hiux | sgi | dec | sun] 


FC 


FORTRAN compiler command name. 


FOPT 


FORTRAN compiler options. 


2. Compilation. 





By executing command make install one executable, i.e. the interface program(grc4f ), 
is generated at BINDIR. Furthermore three libraries, i.e. BASES/SPRING, CHANEL and kine- 
matics utility library, are generated in LIBDIR. 
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3. Install default Makefiles. 

By executing srcO/lgen. sh command, all Makefiles for 76 processes will be generated 
according to the environment where grc4f has been installed. 

The sample control data files will be found at the directory sample. 

TEST RUN OUTPUT 

control data 

process = eNEuD 
end 

Followings are the output files from BASES and SPRING. Only one histogram, the energy 
distribution of the particle 1, is shown since the whole output is too lengthy to be included here. 
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444 ########- 

44444 ### ##- 

###### ##--## ##### 444--44 ########## 

-####---######--##########--####---### 444 44 ### 

-### ## ### ### ## 44444444444444 ### 

— ########## ### ### 44444444444444 ### 

-#### ### ######---### 44 ### 

-############# — ####### ######## 44 ######### 

-## ### 

-############ 

no n r-c tmc mc 



Version 1.0 (2) 
Last date of change : 2 Feb 1996 

Minami*Tateya Collaboration, KEK, JAPAN 
E*mail : grc4f @minami . kek . jp 

Copyright Minami*Tateya Collaboration 

J.Fujimoto, et.al. LEP-II Physics, WW-Generator , CERN 1996. 

Process : e + e- --> (3) electron (4)nu-e-bar (5)u (6)cl-bar 



Canonical CUT: 
Mass (Width) 

W-boson : 
u-quark : 
c-quark : 
t-quark : 



80.230 ( 
.005 
1.300 
170.000 



Z-boson : 
d-quark : 
s-quark : 
b-quark : 



91.189 ( 

.010 

.200 
4 .300 



alpha : 1/128 . 070 
Experimental Cuts 
Angle Cuts 



Particle 


3 




985 


> 




985 


Particle 


4 


: 


000 


> 


: 


000 


Particle 


5 


-l 


000 


> 


: 


000 


Particle 


6 


-l 


000 


> 


: 


000 



Energy CUT 



Particle 3 


: 


000 


> 


190 


000 


Particle 4 




000 


> 


190 


000 


Particle 5 


3 


000 


> 


190 


000 


Particle 6 


3 


000 


> 


190 


000 


Invariant 


mass CUT 












Inv . 


Mass 3-4 




001 


> 


189 


985 


Inv . 


Mass 5-6 


5 


000 


> 


189 


999 


Inv . 


Mass 3-5 




006 




189 


990 


Inv . 


Mass 4-6 




010 




189 


994 


Inv . 


Mass 3-6 




011 




189 


995 


Inv . 


Mass 4-5 




005 


> 


189 


989 



OPTIONS : 

Calculation: TREE 

Anomalous Coupling: NO 

Coulomb Correction: NO 

QCD Correction: NO 

Gluon graph: NO 



only for hadronic process) 



Date: 96/ 2/ 5 17: 
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BASES Version 5.1 
coded by S.Kawabata KEK, March 1994 



<< Parameters for BASES >> 
(1) Dimensions of integration etc. 



# 


of 


dimensions : 


Ndim = 


8 ( 50 at max.) 


# 


of 


Wilds : 


Nwild 


8 ( 15 at max.) 


# 


of 


sample points 


Ncall 


39366 (real) 40000 (given) 


# 


of 


subregions 


Ng 


48 / variable 


# 


of 


regions 


Nregion = 


3 / variable 


# 


of 


Hypercubes 


Ncube 


6561 



(2) About the integration variables 



i 


XL (i) 




xu(i) 


IG(i) 


Wild 


— +- 

1 


.000000E+00 


: 


000000E+00 


1 


yes 


2 


.000000E+00 


: 


000000E+00 


1 


yes 


3 


.000000E+00 


: 


000000E+00 


: 


yes 


4 


.000000E+00 


: 


000000E+00 


: 


yes 


5 


.000000E+00 


: 


000000E+00 


: 


yes 


6 


.000000E+00 


: 


000000E+00 


: 


yes 


7 


.000000E+00 


: 


000000E+00 


l 


yes 


8 


.000000E+00 


: 


000000E+00 


l 


yes 


— +- 













(3) Parameters for the grid optimization step 
Max.* of iterations: ITMX1 = 7 
Expected accuracy : Accl = .1000 % 

(4) Parameters for the integration step 
Max.* of iterations: ITMX2 = 15 
Expected accuracy : Acc2 = .0500 % 



<< Computing Time Information >> 



For BASES 




H 


M 




Sec 


Overhead 













. 04 


Grid Optim. 


Step 





54 





.53 


Integration 


Step 


2 


6 


30 


.43 


Go time for 


all 


3 





31 


.00 



(2 ) Expected event generation time 

Expected time for 1000 events : 18.36 Sec 



Date: 96/ 2/ 5 17:27 
Convergency Behavior for the Grid Optimization Step 



<- Result of each iteration -> <- Cumulative Result -> < CPU time > 

IT Eff R_Neg Estimate Acc % Estimate (+- Error ) order Acc % ( H : M: Sec ) 



: 


18 


.00 


6 


777E- 


01 


2.071 


6 


777356 (+- 


.1403891E- 


01 


2.071 





: 


52 


38 


2 


84 


.00 


6 


928E- 


01 


. 648 


6 


913652 (+- 


. 042747) E- 


CL 


. 618 





9 


52 


98 


3 


86 


.00 


6 


947E- 


01 


.533 


6 


932765 ( + - 


.027989)E- 


01 


.404 





18 


4 


99 


4 


92 


.00 


6 


901E- 


01 


.444 


6 


918370 (+- 


.020668)E- 


01 


.299 





26 


55 


71 


5 


94 


.00 


6 


896E- 


01 


.413 


6 


910548 (+- 


.016720)E- 


01 


.242 





35 


54 


74 


6 


95 


.00 


6 


860E- 


01 


.513 


6 


901197 (+- 


.015105)E- 


o: 


.219 





44 


59 


96 


7 


94 


.00 


6 


981E- 


01 


.783 


6 


906845 (+- 


.014559)E- 


01 


.211 





54 





53 



Date: 96/ 2/ 5 17:27 
Convergency Behavior for the Integration Step 



<- Result of each iteration -> <- Cumulative Result -> < CPU time > 

IT Eff R_Neg Estimate Acc % Estimate (+- Error ) order Acc % ( H : M: Sec ) 



1 88 


.00 


6 


870E- 


01 


.723 


6 


870441 (+- 


.049666)E- 


01 


.723 


: 


2 


24 


93 


2 88 


.00 


6 


918E- 


01 


.723 


6 


894173 (+- 


. 035247) E- 


01 


. 511 


: 


LC 


51 


48 


3 88 


.00 


6 


878E- 


01 


. 710 


6 


888702 (+- 


. 028577) E- 


01 


.415 


: 


19 


17 


98 


4 88 


.00 


6 


897E- 


01 


. 937 


6 


890082 (+- 


.026136)E- 


01 


.379 


: 


27 


44 


75 


5 88 


.00 


6 


868E- 


01 


.669 


6 


884589 (+- 


.022722)E- 


01 


.330 


: 


36 


10 


48 


6 88 


.00 


6 


891E- 


01 


. 695 


6 


885740 (+- 


.020527)E- 


01 


.298 


: 


44 


36 


85 


7 88 


.00 


6 


902E- 


01 


.699 


6 


888267 (+- 


.018890)E- 


o: 


.274 


: 


53 


2 


07 


8 88 


.00 


6 


919E- 


01 


. 708 


6 


892296 (+- 


. 017625) E- 


01 


.256 


2 


: 


28 


62 


9 88 


.00 


6 


900E- 


01 


.709 


6 


893157 (+- 


.016582)E- 


01 


.241 


2 


9 


54 


80 


10 88 


.00 


6 


821E- 


01 


.667 


6 


884725 (+- 


.01558DE- 


01 


.226 


2 


18 


19 


24 


11 88 


.00 


6 


898E- 


01 


. 657 


6 


886173 (+- 


.014734)E- 


01 


.214 


2 


26 


45 


13 


12 88 


.00 


7 


002E- 


01 


. 759 


6 


894471 (+- 


. 014198) E- 


01 


.206 


2 


35 




51 


13 88 


.00 


7 


016E- 


01 


.698 


6 


903859 (+- 


.013637)E- 


01 


.198 


2 


43 


38 


96 


14 88 


.00 


6 


883E- 


01 


. 672 


6 


902177 (+- 


.013080)E- 


01 


.190 


2 


52 


4 


76 


15 88 


.00 


6 


851E- 


01 


.694 


6 


898595 (+- 


.01261DE- 


01 


.183 


3 





30 


96 



********* Contents of the histogram Header ********* 

(1) Actual Buffer size = 9269 Words 

(2 ) Contents of Histograms 

Max . No . of Histograms = 50 
Number of Histograms = 12 



:d 




X_min 




X_max 


X_bin 


Hash 


Hst# 


: 




0000E+00 


: 


9000E+02 


50 


2 


2 


i 


14 




0000E+00 


: 


9000E+02 


50 


2 


2 




2 




0000E+00 


: 


9000E+02 


50 


3 


2 


2 


15 




0000E+00 


: 


9000E+02 


50 


3 


2 


12 


3 




0000E+00 


: 


9000E+02 


50 


4 


: 


3 


4 




0000E+00 


: 


9000E+02 


50 


5 


: 


4 


5 


: 


0000E+00 


: 


0000E+00 


50 


6 


: 


5 


6 


: 


0000E+00 


: 


0000E+00 


50 


7 


: 


6 


7 


: 


0000E+00 


: 


0000E+00 


50 


8 


: 


7 


8 


: 


0000E+00 


: 


0000E+00 


50 


9 


: 


8 


9 




0000E+00 


: 


9000E+02 


50 


10 


: 


9 


LC 




0000E+00 


: 


9000E+02 


50 


1 1 


: 


10 



(3) Contents of Scatter Plots 

Max. No. of Scat_Plots - 50 
Number of Scat_Plots = 



Histogram (ID = 1 ) for Energy of Particle 3 

Linear Scale indicated by "*" 
x d(Sigma)/dx .0E+00 6.2E-03 1.2E-02 1.9E-02 

I E 2 I .000 E 01 I 

I .000 I 2.863+- .960 E -41*000000000000000000000000000 I 

I .038 I 3.338+- .512 E -41*0000000000000000000000000000 I 

I .076 I 3.740+- .532 E -41*0000000000000000000000000000 I 

I .114 I 3.071+- .225 E -41*000000000000000000000000000 I 

I .152 I 3.381+- .147 E -41*0000000000000000000000000000 I 

I .190 I 1.644+- .031 E -31****00000000000000000000000000000000 I 

I .228 I 6.506+- .061 E -31**************00000000000000000000000000000 I 

I .266 I 7.778+- .066 E -31****************0000000000000000000000000000 I 

I .304 I 8.745+- .069 E -31******************00000000000000000000000000 I 

I .342 I 9.791+- .075 E -31********************0000000000000000000000000 I 

I .380 I 1.075+- .008 E -21**********************00000000000000000000000 I 

I .418 I 1.194+- .008 E -21************************0000000000000000000000 I 

I .456 I 1.300+- .009 E -21**************************00000000000000000000 I 

I .494 I 1.441+- .009 E -21*****************************000000000000000000 I 

I .532 I 1.556+- .009 E -21********************************000000000000000 I 

I .570 I 1.707+- .010 E -2I**********************************»0000000000000 I 

I .608 I 1.850+- .011 E -21**************************************0000000000 I 

I .646 I 1.931+- .011 E -21***************************************000000000 I 

I .684 I 1.884+- .011 E -21**************************************0000000000 I 

I .722 I 5.370+- .068 E -31***********0000000000000000000000000000000 I 

I .760 I 4.849+- .223 E -41*00000000000000000000000000000 I 

I .798 I 1.614+- .114 E -41*000000000000000000000000 I 

I .836 I 6.826+- .689 E -51*00000000000000000000 I 

I .874 I 2.196+- .297 E -51*00000000000000 I 

I .912 I 6.076+-2.383 E -61*00000000 I 

I .950 I .000+- .000 E 01 I 

I .988 I .000+- .000 E 01 I 

I 1.026 I .000+- .000 E 01 I 

I 1.064 I .000+- .000 E 01 I 

I 1.102 I .000+- .000 E 01 I 

I 1.140 I .000+- .000 E 01 I 

I 1.178 I .000+- .000 E 01 I 

I 1.216 I .000+- .000 E 01 I 

I 1.254 I .000+- .000 E 01 I 

I 1.292 I .000+- .000 E 01 I 

I 1.330 I .000+- .000 E 01 I 

I 1.368 I .000+- .000 E 01 I 

I 1.406 I .000+- .000 E 01 I 

I 1.444 I .000+- .000 E 01 I 

I 1.482 I .000+- .000 E 01 I 

I 1.520 I .000+- .000 E 01 I 

I 1.558 I .000+- .000 E 01 I 

I 1.596 I .000+- .000 E 01 I 

I 1.634 I .000+- .000 E 01 I 

I 1.672 I .000+- .000 E 01 I 

I 1.710 I .000+- .000 E 01 I 

I 1.748 I .000+- .000 E 01 I 

I 1.786 I .000+- .000 E 01 I 

I 1.824 I .000+- .000 E 01 I 

I 1.862 I .000+- .000 E 01 I 

I E 2 I .000 E 01 I 



x d(Sigma)/dx 1.0E-06 1.0E-05 1.0E-04 1.0E-03 

Logarithmic Scale indicated by "0" 



444 ########- 

44444 ### ##- 

###### ##--## ##### 444--44 ########## 

-####---######--##########--####---### 444 44 ### 

-### ## ### ### ## 44444444444444 ### 

— ########## ### ### 44444444444444 ### 

-#### ### ######---### 44 ### 

-############# — ####### ######## 44 ######### 

-## ### 

-############ 

no n r-c tmc tttc 



Version 1.0 (2) 
Last date of change : 2 Feb 1996 

Minami*Tateya Collaboration, KEK, JAPAN 
E*mail : grc4f @minami . kek . jp 

Copyright Minami*Tateya Collaboration 

J.Fujimoto, et.al. LEP-II Physics, WW-Generator , CERN 1996. 

Process : e+e- --> (3) electron (4)nu-e-bar (5)u (6)d-bar 



Canonical CUT: 
Mass (Width) 

W-boson : 
u-quark : 
c-quark : 
t-quark : 



80.230 ( 
.005 
1.300 
170.000 



Z-boson : 
d-quark : 
s-quark : 
b-quark : 



91.189 ( 

.010 

.200 
4 .300 



alpha : 1/128 . 070 
Experimental Cuts 
Angle Cuts 



Particle 


3 




. 985 


> 




. 985 


Particle 


4 


: . 


.000 


> 


: 


.000 


Particle 


5 


-l. 


.000 


> 


: 


.000 


Particle 


6 


-l. 


.000 


> 


: 


.000 



Energy CUT 



Particle 3 


: 


.000 


> 


190. 


.000 


Particle 4 




.000 


> 


190. 


.000 


Particle 5 


3 


.000 


> 


190. 


.000 


Particle 6 


3 


.000 


> 


190. 


.000 


Invariant 


mass CUT 












Inv . 


Mass 3-4 




.001 


> 


189. 


. 985 


Inv . 


Mass 5-6 


5 


.000 


> 


189. 


.999 


Inv . 


Mass 3-5 




.006 




189. 


. 990 


Inv . 


Mass 4-6 




. 010 




189 . 


.994 


Inv . 


Mass 3-6 




.011 




189. 


. 995 


Inv . 


Mass 4-5 




.005 


> 


189. 


. 989 



OPTIONS : 

Calculation : TREE 

Anomalous Coupling : NO 

Coulomb Correction : NO 

QCD Correction: NO 

Gluon graph : NO 



only for hadronic process ) 
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SPRING Version 5.1 
coded by S.Kawabata KEK, March 1994 



Number of generated events = 


10000 




Generation efficiency = 


16. 


.538 


Percent 


Computing time for generation = 


836. 


.960 


Seconds 


for Overhead = 


1 . 


. 050 


Seconds 


for Others = 




.420 


Seconds 


GO time for event generation 


838. 


.430 


Seconds 


Max. number of trials MXTRY 




50 


per event 


Number of miss-generation = 




56 


times 



************* Number of trials to get an event ************* 
Total = 1005 6 events " * " : No . of events in Linear scale . 





X 


-+■ 


Lg (dN/dx) 


dN/dx . OE+00 6 . 3E+02 1 . 3E+03 1 . 9E+03 




+ ■ 
: 




1 C 


: 




994E 


3 Z 


1 9 94 I **************************************** 000000 0000 


: 




200 


: 


: * 


575E 


3 Z 


1575I ********************************OOOOOOOOOOOOOOOOOI 


: 




300 


: 


: * 
- * 


249E 


3 Z 


12 4 9l *************************ooOOOOOOOOOOOOOOOOOOOO 


I 


: 




400 


: 




3 60E 


2 : 


9 36I *******************ooooooooooooooooooooooooooo 


I 


: 




500 


: 


7 ■ 


3 90E 


2 : 


73 91***************00000000000000000000000000000 


I 


: 




600 


: 


5 , 


77 ? E 


2 : 


577I ************ OOOOOOOOOOOOOOOOOOOOOOOOO OOOOO 


I 


: 






: 


4 ' 


E 


2 : 


4481*********00 OOOOOOOOOOOOOOOOOOOOOOOOO OOOOO 


I 


: 




«nn 


: 


4 ' 




2 : 


4161********* OOOOOOOOOOOOOOOOOOOOOOOOOO OOOOO 


I 


: 






: 




0E 




3021****** * 0000000000000000000000000000000 


I 


: 




000 


: 


2 , 


650E 


. : 
2 : 


2651****** OOOOOOOOOOOOOOOOOOOOOOOOOO OOOOO 


I 


: 


: ■ 


Z 


: 


2 . 


1 90E 


2 : 


2191*****0000000000000000000000000000000 


I 


: 


: ■ 


200 


: 




910E 


2 : 


1911****0000 000000000000000000000000000 


I 


: 


: ■ 
: ■ 


300 


: 


: ■ 
: ■ 


530E 


2 : 


1 5 3 1 * * * * 000000000000000000000000000000 


I 


: 


: ■ 


400 


: 


: * 


330E 


2 : 


1331***000000000000000000000000000000 


I 


: 


: ■ 


500 


: 


: ■ 


000E 


2 : 


1001***0000000000000000000000000000 


I 


: 




600 


: 


- ■ 


130E 


2 : 


1 1 3 1 * * *00000000000000000000000000000 


I 


: 


: ■ 


700 


: 


7 ■ 


800E 




781**00 00000000000000000000 OOOOO 


I 


: 


: ■ 


800 


: 


6 , 


200E 


:: 


621 * * OOOOOOOOOOOOOOOOOOOOOOOOOO 


I 


: 


: ■ 


900 


: 


5 . 


900E 


:: 


5 91**0000000000000000000000000 


I 


: 


: ■ 

2 . 


000 


: 


4 ' 


700E 


:: 


4 71*0000000000000000000000000 


I 


: 


2 . 


100 


: 


4 ' 


200E 


:: 


4 2 1*000000000000000000000000 


I 


: 


2 . 


200 


: 


3 . 


300E 


:: 


331*000 0000000000000000 000 


I 


: 


2 . 


300 


: 


2 . 


700E 


:: 


271*00 000000000000000 000 


I 


: 


2 . 


400 


: 


3 , 


500E 


:: 


351*00000000000000000000000 


I 


: 


2 . 


500 


: 


2 . 


700E 


:: 


271*000000000000000000000 


I 


: 


2 . 


600 


: 


2 , 


000E 


:: 


201*0000000000000000000 


I 


: 




700 


: 


2 . 


300E 


:: 


2 3 1 * 000000000000000 OOOOO 


I 


: 


9 


800 


: 


" ' 


500E 


:: 
_ " 


151*0 00000000000 0000 


I 


: 


2 . 


900 


: 


9 . 


000E 


_ 


91*00000000000000 


I 


: 


3 , 


000 


: 


2 , 


000E 




201*0000000000000000000 


I 


: 


3 , 


100 


: 


9 . 


000E 


c _ 


91*00000000000000 


I 


: 


3 , 


200 


: 


7 ' 


000E 


z 


71*000000000000 


I 


: 


3 , 


300 


: 


7 " 


oooe 


z 


71*000000000000 


I 


: 


3 . 


400 


: 


" ' 


100E 




111*000000000000000 


I 


: 


3 , 


500 


: 


6 . 


oooe 


c _ 


61*00000000000 


I 


: 


3 , 


600 


: 


6 , 


oooe 


z 


61*00000000000 


I 


: 


3 , 


700 


: 


5 . 


OOOE 


c z 


51*0000000000 


I 


: 


3 . 


800 


: 


B . 


OOOE 


°" 


8 1 * 0000000000000 


I 


: 


3 , 


900 


: 


4 ' 


000E 


c _ 


41*000000000 


I 


: 


4 ' 


000 


: 


2 . 


000E 


z 


21*0000 


I 


: 


4 ' 


100 


: 


3 . 


000E 


c z 


31*0000000 


I 


: 


4 ' 


200 


: 


4 ' 


000E 
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Lg (dN/dx) 


dN/dx 1. OE+00 l.OE+01 1.0E+02 
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"0" : No . of Events in Log . scale . 



********* contents of the histogram Header 

(1) Actual Buffer size = 9269 Words 

(2) Contents of Histograms 

Max . No . of Histograms = 50 
Number of Histograms = 12 



ID 




Xmin 






X_max 


X_bin 


Hash 


Hst# 


: 
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: . 


9000E+02 


50 


2 


2 


l 


14 
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Scatter Plots 










Max 


. No. of 


Scat_ 


_Plots = 
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Number of 


Scat. 


_Plots - 












Original Histogram (ID = 1 ) for Energy of Particle 3 
Total - 10000 events "*" : Orig. Dist. in Log Scale. 



d(Sig/dx) dN/dx 1.0E-06 1.0E-05 1.0E-04 1.0E-03 



I E 2 I .000E 01 01 I 

I .000 I 2.863E -41 4 I*****************o*0*>***** I 

I .038 I 3.338E -41 14 I ...*.***.*.*.*.*.*.***.*.< >* I 

I .076 I 3.740E -41 9 i**********************oo>*** j 

I .114 I 3.071E -41 I4i*************************< > I 

I .152 I 3.381E -41 14 I*************************<0>* I 

I .190 I 1.644E -31 841**********************************0 I 

I .228 I 6.506E -31 380 I******************************************o I 

I .266 I 7.778E -31 38 0I******************************************o* I 

I .304 I 8.745E -31 54 71********************************************0 I 

I .342 I 9.791E -31 5 08I*******************************************o I 

I .380 I 1.075E -21 60 6I********************************************o I 

I .418 I 1.194E -21 6891*********************************************0 I 

I .456 I 1.300E -21 6 971*********************************************0 I 

I .494 I 1.441E -21 784 1*********************************************0 I 

I .532 I 1.556E -21 85 61**********************************************0 I 

I .570 I 1.707E -21 9641**********************************************0 I 

I .608 I 1.850E -21 10 61***********************************************0 I 

I .646 I 1.931E -21 10371***********************************************0 I 

I .684 I 1.884E -21 104 81***********************************************0 I 

I .722 I 5.370E -31 3061*****************************************0 I 

I .760 I 4.849E -41 341******************************0 1 

I .798 I 1.614E -41 i 3 i************************< > I 

I .836 I 6.826E -51 51******************0*0 > I 

I .874 I 2.196E -51 n************* * > j 

I .912 I 6.076E -61 I********* I 

I .950 I .000E 01 01 I 

I .988 I .000E 01 01 I 

I 1.026 I .000E 01 01 I 

I 1.064 I .000E 01 01 I 

I 1.102 I .000E 01 01 I 

I 1.140 I .000E 01 01 I 

I 1.178 I .000E 01 01 I 

I 1.216 I .000E 01 01 I 

I 1.254 I .000E 01 01 I 

I 1.292 I .000E 01 01 I 

I 1.330 I .000E 01 01 I 

I 1.368 I .000E 01 01 I 

I 1.406 I .000E 01 01 I 

I 1.444 I .000E 01 01 I 

I 1.482 I .000E 01 01 I 

I 1.520 I .OOOE 01 01 I 

I 1.558 I .OOOE 01 01 I 

I 1.596 I .OOOE 01 01 I 

I 1.634 I .OOOE 01 01 I 

I 1.672 I .OOOE 01 01 I 

I 1.710 I .OOOE 01 01 I 

I 1.748 I .OOOE 01 01 I 

I 1.786 I .OOOE 01 01 I 

I 1.824 I .OOOE 01 01 I 

I 1.862 I .OOOE 01 01 I 

I E 2 I .OOOE 01 01 I 



x d(Sig/dx) dN/dx "0" : Generated Events.! Arbitrary unit in Log ) 



